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Abstract 

We generate new hierarchy of many-parameter family of maps of the interval [0, 1] 
with an invariant measure, by composition of the chaotic maps of reference Jl|]. Using 
the measure, we calculate Kolmogorov-Sinai entropy, or equivalently Lyapunov char- 
acteristic exponent, of these maps analytically, where the results thus obtained have 
been approved with numerical simulation. In contrary to the usual one-dimensional 
maps and similar to the maps of reference p]] , these maps do not possess period dou- 
bling or period-n-tupling cascade bifurcation to chaos, but they have single fixed point 
attractor at certain region of parameters values, where they bifurcate directly to chaos 
without having period-n-tupling scenario exactly at these values of parameter whose 
Lyapunov characteristic exponent begins to be positive. 

Keywords: Chaos, Invariant measure, Entropy, Lyapunov characteristic ex- 
ponent, Ergodic dynamical systems . 
PACs numbers:05.45.Ra, 05.45.Jn, 05.45.Tp 
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In the past twenty years dynamical systems, particularly one dimensional iterative maps 
have attracted much attention and have become an important area of research activity. One 
of the landmarks in it was introduction of the concept of Sinai-Ruelle-Bowen (SRB) measure 
or natural invariant measure 0, H|. This is, roughly speaking, a measure that is supported 
on an attractor and also describe the statistics of the long time behavior of the orbits for 
almost every initial condition in the corresponding basin of attractor . This measure can be 
obtained by computing the fixed density of the so called Frobenius-Perron operator which 
can be viewed as a differential-integral operator, hence, exact determination of invariant 
measure of dynamical systems is rather a nontrivial task, such that invariant measure of 
few dynamical systems such as one-parameter family one-dimensional piecewise linear maps 
[[J. |5[ including Baker and tent maps or unimodal maps such as logistic map for certain values 
of its parameter, can be derived analytically. In most cases only numerical algorithms, as 
an example Ulam's method |6|, |7], || are used for computation of fixed densities of Frobenius- 
Perron operator. 

Authors in reference [JJ have given hierarchy of one parameter family of nonlinear maps 
of interval [0, 1] with an invariant measure. Here in this paper we generate new hierar- 
chy of many-parameter family of maps of the interval [0, 1] with an invariant measure, by 
composition of the chaotic maps of reference Jl|. These maps are also defined as ratio of 
polynomials, where we have derived analytically their invariant measure for arbitrary values 
of the parameters. Using this measure, we have calculated analytically, Kolmogorov-Sinai 
entropy or equivalently positive Lyapunov characteristic exponent of these maps, where the 
numerical simulation approve the analytic calculation. Also it is shown that just like the 
maps of reference they possess very peculiar property, that is, contrary to the usual, 
these maps do not possess period doubling or period-n-tupling cascade bifurcation to chaos, 
but instead they have single fixed point attractor at certain region of parameters values, 
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where they bifurcate directly to chaos without having period-n-tupling scenario exactly at 
these values of parameter whose Lyapunov characteristic exponent begins to be positive. 
The paper is organized as follows: In section II we introduce new hierarchy of family of 
many-parameters maps by composition of the chaotic maps of refernce 0. In Section III we 
show that the proposed anzats for the invariant measure of these maps are eigenfuntion of 
Ferobenios- Perron operator with largest eigenvalue 1. Then in section IV using this measure 
we calculate Kolmogorov- Sinai entropy of these maps for an arbitrary value of parameters. 
In section V we compare analytic calculation with the numerical simulation. Paper ends 
with a brief conclusion. 



2 Many- Parameter Families of Chaotic Maps 

Let us consider the one-parameter families of chaotic maps of the interval [0, 1] given in 
reference defined as the ratio of polynomials of degree iV: 

a 2 (l + (-l) N 2 F 1 (-N,N,lx) 



a 2 + 1) + (a 2 - 1)(—1) N 2 F 1 (—N, N, \, x) 
a 2 {T N {^)f 



2 ' 

.,2(rn t rZ\\2 



(2-1) 



2 F 1 (-N, N,-,x) = (-1)^008 (2A^arccos v / ^) = (-l) N T 2N (y/x) 



l + (a 2 -l)(T N (^) 2 ) 
where iV is an integer greater than one. Also 

1 

2' 

is hypergeometric polynomials of degree N and Tn(x)(Un{x)) are Chebyshev polynomials 

of type I (type 11)0, respectively. Obviously these map the unit interval [0,1] into itself. 

$tv(x, a) is (N-l)-model map, that is it has (N — 1) critical points in unit interval [0, l],(see 

Figure 4) since its derivative is proportional to derivative of hypergeometric polynomial 

2 Fi{— N, N, |, x) which is itself a hypergeometric polynomial of degree (N — 1), hence it has 

(N — 1) real roots in unit interval [0, 1]. Defining Shwarzian derivative ||10|| SQn(x,oi) as: 

Qfd. { )) m^*) 3/ n(x,cQ \ 2 ( *» N (x,a) Y l( ^(x,a) \ 2 
{ N{ ' )] ®' N {x,a) 2\& N (x t a)J \& N {x,a)) 2\V N (x,a)J ' 
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with a prime denoting a single differentiation with respect to variable x, one can show that [[J: 

S ($ N (x, a)) = S (zFti-N, N, i x)) < 0. 



Therefore, the maps &n(x) have at most N + 1 attracting periodic orbits[ID]. As it is shown 



in reference [[IJ, these maps have only single period one stable fixed points. 
Using the above hierarchy of family one-parameter of maps we can generate new hierarchy 
of family many-parameters chaotic maps with an invariant measure simply from the compo- 
sition of these maps. Hence considering the functions ^^(x, a k ), k — 1, 2, • • • ,n we denote 

Nl,N 2 ,-,N n 



their composition by: $^'^ 2 ' ".'.' a N n ( x ) which can be written in terms of them in the following 



form: 



$ Nl ($ Ni (- ■ ■ {$N n (x, «„),«(„_!)) • • • , a 2 ),a 1 ) (2-2) 

Since these maps consist of the composition of the (N k — l)-modals (k = 1, 2, • • • , n) maps 
with negative Shwarzian derivative, therefore, they are (N\N2---N n — 1)— modals map 
and their Shwarzian derivative is negative too fl0| . Therefore these maps have at most 
N1N2 ■ ■ ■ N n + 1 attracting periodic orbits [|IU||. As we will show below in this section, these 
maps have only a single period one stable fixed points. Since, denoting m-composition of 
these functions by it is straightforward to show that the derivative of $ ( - m ) at its possible 
m x n periodic points of an m-cycle: x^ k+1 = $^(3^*., a k ), x ljfi+1 = § Nn (x n ^, a N ), ll = 
1,2 ■■■,m and k — 1, 2, •■• ,n and x^x — ^N n {x m ,m a n) is 

1 ^ (m) \== noi 1 — («2 + (i - i' ( 2 - 3 ) 

ax ^=1 fc=i a ^ 

since for G [0, 1] we have: 

min(al + (1 — a 2 k x^ k )) = mm(l, a^), 

therefore, 

min I .. , ± 
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Hence the above expression is definitely greater than one for n k=i jtt~ < ITfc=i a k < Uk=i Nk, 
that is, these maps do not have any kind of m-cycle or periodic orbits in the region of the 
parameters space defined by Y\ k=1 < Uk=i a k < Uk=i Nk, actually they are ergodic in this 
region of the parameters space. From (2-3) it follows that | £$ {m) | at mxn periodic points 
of the m-cycle belonging to interval [0,1], varies between Y\ k=1 (N k a k ) m and YYk=\ {—) f° r 
nLi«* < IILi ]fe and between nLi(f) m and Y[ n k=l {N k a k ) m for nLi > \Tk =1 N k , 
respectively. 

From the definition of these maps, we see that definitely x — 1 and x = (in special case of 
odd integer values of N — 1, N 2 , • • • , N n ) belong to one of the m-cycles. 
For Uk=i a k < Uk=i j^(Ilfc=i a k > Uk=i N k), the formula (2 - 3) implies that for those cases 
in which x = l(x = 0) belongs to one of m-cycles we will have | -jk^™^ \< 1, hence the curve 
of $( m ) starts at x — l(x — 0) beneath the bisector and then crosses it at the previous (next) 
periodic point with slope greater than one (see Fig. 1), since the formula (2 — 3) implies 
that the slope of fixed points increases with the decreasing (increasing) of | x^ k |, therefore 
at all periodic points of n-cycles except for x — l(x — 0) the slope is greater than one that 
is they are unstable, this is possible only if x — l(x — 0) is the only period one fixed point 
of these maps. 

Hence all m-cycles except for possible period one fixed points x — 1 and x = are unstable. 
Actually, the fixed point x = is the stable fixed point of these maps in the regions of the 
parameters spaces defined by a k > 0,k = 1, 2, • • • , n and H k= 1 a k < H k=1 ^- only for odd 
integer values of N 1: N 2 , ■ ■ ■ , N n , however, if one of the integers N k , k = 1.2, • • • , n happens 
to be even, then the x = will not be a stable fixed point anymore. But, the fixed point 
x — 1 is stable fixed point of these maps in the regions of the parameters spaces defined by 
]lfc=i otk > Uk=i Nk and a k < 00, k — 1, 2, • • • , n for all integer values of N u N 2 , ■ ■ ■ , N n . 
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As an example we give below some of these maps: 



ai 2 Ux(x-1) + (2x-l) 2 a 2 2 Y 

02 2 Q2 (^) = " " " " t- (2-4) 

«i 2 (Ax(x- 1) + (2x- lfa 2 2 ) +hl 



? 2,3 



? 3,2 



iai,«2 , 



ai 2 ((x - 1) (4 x - if + x (4 x - 3) 2 a 2 2 ) 2 
«i 2 ((x - 1) (4x - if + x (4 x - 3) 2 « 2 2 ) 2 + /i2 

= « 2 2 ((x-l)(4x-l) 2 + x(4x-3) 2 a 1 2 ) 2 
a 2 2 ((x-l)(4x-l) 2 + x(4x-3) 2 a 1 2 ) 2 + /i3 

ai 2 a 2 2 x (4 x - 3) 2 (3 (x - 1) (4 x - l) 2 + x (4 x - 3) 2 a 2 2 ) 

' 3 ' 3 VXj ~ - (x - l) 3 (4 x - l) 6 + 3 x (3 ai 2 - 2) (4 x - 3) 2 (x - if {Ax- if a 2 2 + hi 

(2-7) 

where : 

hi = -16a 2 2 (2x- l) 2 x(x- 1) 

hi = -4x(x - l)(4x - l) 2 (4x - 3) 2 a 2 2 

h3 = 4x(x - l)(4x - l) 2 (4x - 3) 2 «i 2 

hA = 3 x 2 (-3 + 2 af) {x - 1) (4 x - l) 2 {Ax- 3) 4 a 2 4 + c^x 3 (4 x - 3) 6 a 2 % 

Below we also introduce their conjugate or isomorphic maps which will be very useful in 
derivation of their invariant measure and calculation of their KS-entropy in the next section. 
Conjugacy means that the invertible map h{x) = maps / = [0, 1] into [0, 00) and 
transforms maps $jv fc (x, at) into $Ar fe (x, a.k) defined as: 

$jv fc (x, otk) = h o $jv fc (x, a k ) o /i* -1 ) = — tan 2 (iV fc arctan y/x) 
Hence this transforms the maps ^^^ 2 ['.^ n {x) into ^^^'['.'° l N n ( x ) defined as: 
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—2 tan 2 (A r 1 arctan v^°^j tan 2 (^2 arctan • • • — ^ tan 2 (iV n arctan \/^))) 



at 



Ho 



a? 



tan 2 (A^! arctan 



\ 



«2 



tan 2 (A^ 2 arctan 



0'" 



tan 2 (A^„ arctan v/x) •• •)) (2-8) 



3 INVARIANT MEASURE 

Dynamical systems, even apparently simple dynamical systems as those described by maps 
of an interval, can display a rich variety of different asymptotic behavior. On measure 
theoretical level these types of behavior are described by SRB || or invariant measure 
describing statistically stationary states of the system. The probability measure //on [0, 1] 
is called an SRB or invariant measure of the maps $^'m 2 '".'. , .'jv ( x ) gi ven i n (2 — 2), if it is 
^iVi'iV2 (^-invariant and absolutely continuous with respect to Lebesgue measure. For 
deterministic system such as these composed maps, the $^'^'.^ n (:r)-invariance means that 
its invariant measure u^au^.-^n (x) fulfills the following formal Ferbenius-Perron integral 
equation 

This is equivalent to: 

/',,,":•".'• ••'. (y) = fl^W-^n (x)^- , (3-1) 

If $ a l- a 2.-.«« ( v ) 12'™ U,y 

defining the action of standard Ferobenius-Perron operator for the map $(x) over a function 



as: 



dx 

P^ T --, 7 f(y)= £ f(x)- . (3-2) 
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We see that, the invariant measure u*«i.«2.-.«» (x) is given as the eigenstate of the Frobenius- 

r *JV 1 ,JV 2l ...,JV„ 

Perron operator P*ai,a 2 .-.«„ corresponding to the largest eigenvalue 1. 

N 1 ,N 2 ,---,N„ 

As we will prove below the measure u s ai,«2.-.«„ defined as: 

^AT 1 ,iV 2 ,...,JV n 



^ (3-3) 



is the invariant measure of the maps ^^ a N^..^N n ( x ) P rov ided that the parameter (3 is positive 
and fulfills the following relation: 



k R TT\ R f 1 ^ R / 1 ~ X ' X r 7 1 ~ - 1 W" 4 J 

fc = l ^NnKg) J->N n -l\ a n -DjV n _ 2 l a n _!,a„ x J V « 2 ,a 3 , •• / J_ J 

P *?iV„^ r 'jV„_ 1 ,JV„^i r 'jV 2 ,JV 3 ,...,JV„i / 3' 1 



where the polynomials A Nk (x) and B Nk (x) (k — 1, 2, • • • , n) are defined as: 



where [ ] means greatest integer part. Also the functions VN^ip), VnI^nS^) ■> ' ' ' an< ^ 
Vn^n^'-^nS^) are defined in the following form: 



1 a n _iA A r n _ 1 ( -5^ i ) 
a„_i,a„ / \ />/ 'iN n \p> \2 

VNn-uNnKlS) - P\ R ( 1 \ J 



1 Q^-^Afa („ a 3. a 4.---. c «n TTn ) 

a2,a 3 ,-,a„ / N 'N 3 ,N 4 ,---,N n y p > ^2 

r lN2,N3,-,N n \ a) — P\ r> ( 1 \ J 

M J ->N2\„ a 3' a 4' — ,a n 

'lN 3 ,N 4 ,---,N n [ • (3 1 

As we see the above measure is defined only for (3 > hence, from the relations (3 — 4), 
it follows that these maps are ergodic in the region of the parameter space which leads to 
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positive solution of (3. Taking the limits of (3 — > + and (3 — > oo in the relation (3- 
4), respectively, one can show that the ergodic regions are : YYl=i ^ < Uk=i a k < Uk=i Nk 
for odd integer values of Ni, N 2 , ■ ■ ■ , N n and a k > 0, for A; = 1, 2, • • • , nk, Y\!k=i a k < Uk=i Nk 
if one of the integers happens to become even, respectively. Out of these regions they have 
only single stable fixed points. 

In order to prove that measure (3 — 3) satisfies equation (3 — 1), with (3 given by relation 
(3 — 4), it is rather convenient to consider the conjugate map: 

®n u N2,-,n„oci, a 2 , ■ ■ ■ , a n (x), 



with measure ui«i,«2-,«» denoted by h-, related to the denoted by us 

through the following relation: 

Denoting $jVi,Jv 2 ,-,iv„Q!i, «2, • • • , os n (x) by y and inverting it, we get : 

x kl = tan 2 (— arctan Jt/a:? + — ) h = 1,.., TVi. 

2 1 , / ~ t k 2 7l. 



x klM =tan (— arctan J x kl a?i + — ) k 2 = l,..,N 2 . 
l\ 2 l\ 2 



x kl M,-,k n = tan (tt arctan x klM ... >kn _ 1 a 2 n + — ) h = 1, .., JVi. 



Then, taking derivative of x klt k 2 ,.-,k n with respect to y, we obtain: 



dx k i t k 2 ,-,k„ _ j 




1 x ki,k2,---,k„ 


dy 







(1 + ^fcl,fc 2 ,-,fcn)(l + X k 2 ,k 3 ,-,kn) •••(! + ^fc„-l,fcn)(l + Xfc n ) . 

(1 + a^ fe2)jt3i ... jfen )(l + an_i^ 3 ,fc4,-,fcn) •••(! + ai^„_ 1 ,fe„)(l + & 2 2 x kn )(l + a?y) ' 
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In derivation of the formula we have used the chain rule properties of the derivative of 
composite functions. 

Substituting the above result in equation (3 — 1), we have: 

k=l k ki &2 k n 



n )(1 + %k2,kz,---,k (1 + 

)(1 + X kn ) 

(1 + alx k2M ,... >kn ){l + a£_iXfc 3 ,fc4,-,fcJ ■-•(! + «i^n-i,fe»)( 1 + ai^fcj 
Now, considering the following anzatz for the invariant measure jx$(y): 



then the above equation reduces to: 

l + Py [ t}i N k 

x y y ...y ( ^ ± x ^M,-,k n )H + x k2M ,., kn ) ■••(! + x kn _ ukn )(i + x kn ) \ 

Now, using the formula (3-5) of reference [l]] we obtain: 

j2 1 + x kl ,k 2 ,-,k n _ a nA Nn (f) 1 + a>lx hlMt „., fcB _ 1 

I, i 1 



a, 



AT„ fcn=1 1 + (3x kl ,k 2 ,-,k n B Nn {^) 1 + VNli^kiM^kn-i 

a n-l a n ^ ( 1 + ^fci.fca,-,*^-! ) (1 + ,k2,-,k n ) 



A^n-lA^n fcn _ 1 = lfen = 1 (l + «n-l^fci,fca I .» ) *n_i)(l + ^fci,fc 3 ,-,fc n 

a n o; n _i Ajv n (|) Ajv ( i. ) i , fl 2 

X 



n Ni «2 

IIj)EE 



fc=i 




(1 + £ fcl , fc2 ,.., fc J(l + ^^fcj ■•■(! + a; fc „_ 1[fc „)(l + XfcJ 
(1 + a£x fa , fc8l ... )fc J( 1 + a 2 n _xX kzMr . tkn ) •■■(! + «i a; fc„-i,fe„)( 1 + ai^fcn), 
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X TTT X 



k=l B N n {^) Bn^^jt^) 

A-Nn-A "n-wn ,1 J A Nl ( a 2 ,a 3 ,.?".,o„ , \ , ) i , 2,. 



X ; — X • X 



Now, inserting the right hand side of last relation in (3 — 5), we get: 

A-Nn-2 ( «n-l.«n fi J ^JVi ( «a.*3."-.*n. f 1 s ) i , 2„ 

X * X 



We see that the above relation will hold true provided that the parameter /? fulfills the re- 
lation (3-4). 



4 KOLMOGROV-SINAI ENTROPY 

Kolomogrov- Sinai entropy (KS) or metric entropy measure how chaotic a dynamical sys- 
tem is and it is proportional to the rate at which information about the state of dynamical 
system is lost in the course of time or iteration. Therefore, it can also be defined as the aver- 
age rate of information loss for a discrete measurable dynamical system # 2 '"?°lv (#) , ^) , 

by introducing a partition a = A c (ni, n 7 ) of the interval [0, 1] into individual laps A,i one 

can define the usual entropy associated with the partition by: 

71(7) 

H(/i, 7) = - m ( AO In m(A c ) , 

i=l 
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where m(A c ) = J ne A i fJ>(x)dx is the invariant measure of Ai. Defining n-th refining 7(71) of 
7: 

7 B = U(nx:#»(*))~ (fc H7), 

and defining an entropy per unit step of refining by : 

K^^nM^) = lim_ 00 (i//(/i,7)), 

if the size of individual laps of "f(N) tends to zero as n increases, then the above entropy is 
known as Kolmogorov-Sinai entropy that is: 

*%$£%( x )) = Ma». ^1, ^2, • • • , ^' a2 '-' a "x), 7 ). 

KS-entropy , which is a quantitative measure of the rate of information loss with the refining, 
may also be written as: 

Kit, = J 1 ^naw i' (4- 1 ) 

which is also a statistical mechanical expression for the Lyapunov characteristic exponent, 
that is, the mean divergence rate of two nearby orbits. The measurable dynamical system 
(<&iVi, N 2 , • • • , N n ? u0t2 ''"' an x), /i) is chaotic for h > and predictive for h — 0. 
In order to calculate the KS-entropy of the maps ^^'^".'.'.^(aO, it is rather convenient to 
consider their conjugate maps given by (2 — 8), since it can be shown that KS-entropy is a 
kind of topological invariant, that is, it is preserved under conjugacy map. Hence we have: 

Using the integral (4 — 1), the KS-entropy of ^^'^".'.'.^ (x) can be written as 
1 f 00 yffidx d 1 

- / 7=7; T ln (l -r (-2 tan 2 (AT ia rctan y/y N2 ,N 3 ,.,N n ))x 
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d 1 d 1 

{— tan 2 (A^ 2 arctan ^yN 3 ,N 4 ,-,N n )) • • • ^~(— tan 2 (7V n arctan v 7 ^)) |) 



or 



where 



y/j3dx d 1 



/■oo y pdx d 1 

/ i /? ^ ln (l 1 (— tan 2 (^ arctan ^yN 2 ,N 3 ,-,N n ))) + 

1 r°° y/Pdx d 1 



Wo y/x(l + Px) dy N ^ N3 ,;N n «i 
/■oo \J~j3dx d 1 

/ i ft S ln d 1 (~a tan 2 !^ arctan y/y N3 ,Ni,.,N n ))) + ••■ + 

n Jo y/x{l + (3x) dy Ns:N4i . :Nn at% 



1 f°° J]3dx . d . 1 o/ a, , x\x 

"/ / — i ft ^ 1 (-2— tan (iV n _i arctan v /^)))+ 

1 f°° y/ddx . .. d , 1 

U v^OTM 111(1 Tx { al 



' ^ff. , Ml -f(^tan 2 (AT n arctan^)))) |) (4-2) 
jo \/x(l + ox) dx at 



VN n = ^-tan 2 (A? n arctan( v /a : )) (4-3) 



a 2 



2/JV n _i,JV„ = ^^tan 2 (iV n _iarctan( v /^))) (4-4) 



yN 2 ,N 3 ,-,N n = \ tan 2 (A^i arctan(^/yjv3 i jv 4 ,-,jv n ))- (4-5) 

. Now, we calculate the integrals appearing above in the expression for the entropy, sepa- 
rately. The last integral in (4-2) can be calculated with the prescription of section IV of the 
reference Q and it reads: 



n Jo y/x(l + px) dx a 2 n " ' I A Nn (±)B Nn (±) 



(4-6) 
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In order to calculate the integral one before last in (4-2), that is: 

~ I" r^^R \ ln « ^-(-^-tan^JV^arctanVi^))) (4-7) 
7r 7o Vx(l + fix) dy Nn o? n _ x 

first we make the following change of variable by inverting the relation(4-3): 

^fc n _! = tan (— — arctan( v ^;a n _ 1 ) + — ) k n - X = 1, ••, A^-i- 

then the integral (4-7) reduces to: 

£ - /. fcn_1 . =4tT7 T V ln (l tan 2 (iV n _i arctan ^Tj) 

where ^fc n _ 1 and aj£ _ (fc n _i = 1, 2, ■ • • , iV n _i) denote the initial and end points of fc-th branch 
of the inversion of function y^ n = (— tan 2 (iV n arctan y/x)), respectively. Now, inserting the 
derivative of Xk n _ 1 with respect to yN n in the above relation and changing the order of sum 
and integration, we get: 

1 r° fZi a n _iyxT~7(l + x fcn _ 1 ) 



vr Jo k ^ =l v n N n - ly /y^{l + <_i2/iv n ) V^TC 1 + P x kn-i) 

xln(| —^—(^^tan 2 (A^ n „i arctan,/^))). 

Using the formula (3-5) of reference ^ reduces to 

1 [°° y/fidy n B Nn {±) a n A Nn {f) d 1 2 

- / — ^ , m + „ 7TY ln (l 3 (~2— tan (^"-i arctan v 7 ^)))- 

Finally, calculating the above integral with the prescription of reference M we obtain: 



In 



'N n _ l (l + V Z+2^f n ) N ^- i 



Similarly, we can calculate the other integrals appearing in the expression for the entropy of 
the composed maps given in (4-2): 

1 f°° y/j3dx , d ,1 2/Ar , ., 

vr Jo y/x(l + Px) dy NkjNk+lj . jNn a l k _ x 
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A (<Xk,<Xk+i,-,an/l\\T3 ( a k ,a k+1 ,-,a„ w n J.tOr/C 1, Z, • • • , U 

A N k _ 1 {VN k ,N k _ u ---,N n Kp))&N k _AVN k ,N k _ u ---,N n ){p) J 

Finally summing the above integral we get the following expression for the entropy of these 
maps: 



(M ■ ■ ■ iV w )(l + vgg^Hl + yggg)^^ " • " (1 + ygggMlj!^ 

(4-£ 



Imposing the relations between the parameters = 1, 2, • • • , n which are consistent with 
the relation (3-4), reduces these maps to other maps of many-parameters family of the maps 
with number of the parameters less than n. Particularly by imposing enough relations we 
can reduce them to one-parameter family of chaotic maps with an arbitrary asymptotic 
behavior as the parameter takes the limiting values. Hence we can construct chaotic maps 
with arbitrary universality class. As an illustration we consider the chaotic map ^^(a?)- 
Using the formula (4-8) we have 

with the following relation among the parameters a±, a 2 and (5: 

e*i(4/3 + a 2 2 (l + (3f) = 4a 2 /3(l + (3) 

which is obtained from the relation(3-4). Now choosing f3 = a 2 , , < v < 2, the above 
relation reduces to: 

4a£ + "(l + a%) 



ai 



al(l + a£) 2 + 4a£ 



and entropy given by (4-8) reads: 

(l + a|) 2 (2«! + a 2 (l + a»)) 2 



M/i,$? 2 2 (x) = ln. 



;i + a 2 -)(4^ + ai(l + ^) 2 ) 
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which has the following asymptotic behavior near a 2 — ► and a 2 — ► oo: 

h([J,, $22 ( X ) ~ a 2 asa 2 ► 

$2,2^) ~ (^)* as «2 ► OO. 

The above asymptotic form indicates that, for an arbitrary value of < v < 2, the maps 
$2 2( a; ) belong to the universality class which is different from the universality class of chaotic 
maps of the reference |IJ or the universality class of pitch fork bifurcating maps. 



Here in this section we try to calculate Lyapunov characteristic exponent of maps <& N ' (x, a), 
N = 1,2,... .,5 in order to investigate these maps numerically. In fact, Lyapunov charac- 
teristic exponent is the characteristic exponent of the rate of average magnificent of the 
neighborhood of an arbitrary point xq and it is denoted by A(xq) which is written as: 



^Ni,Na,-,N n \ x 0) — tWln-^oo | ^ N 1 ,N 2 ,-,N n { X ^ a > ^Ni,N 2 ,~,N n — - ° W Ni,N 2 ,-,N n I 

= lim n ^ £ In | ^ 2 X" ( fc ' } I, (4-1) 



fc=0 



where x fc = ° ^X'-'X ° - ° ■ lt is obvious that A (1 ' 2) (^o) < for 

an attractor, AN U N 2 , ■ ■■ , N n ai ' a2 '-' an (x ) > for a repeller and AN h N 2 ,---, N n ai ' a2 >~' an (x ) = 
for marginal situation. Also the Liapunov number is independent of initial point Xq, pro- 
vided that the motion inside the invariant manifold is ergodic, thus AN\, N 2 , ■ ■ ■ , N n ai ' a2 ''"' an {xo) 
characterizes the invariant manifold of ^jvV^'-'-'iVn as a wn ole. For values of parameter 
ak,k = 1, 2, • ■ ■ , n, such that the map $nV,n 2 \--'^v„ be measurable, Birkohf ergodic theorem 
implies equality of KS-entropy and Liapunov characteristic exponent, that is: 

Kv, = A^;:::X(x , ( 4 - 2 ) 

Comparison of analytically calculated KS-entropy of maps §^^]'' a fi n (x, a) for N x = 2,3 
and N 2 = 2,3, (see Figures 5, 6 and 7 ) with the corresponding Lyapunov characteristic 
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exponent obtained by simulation, indicate that in chaotic region, these maps are ergodic as 
predicted by Birkohf ergodic theorem. In non chaotic region of the parameter, Lyapunov 
characteristic exponent is negative, since in this region we have only stable period one fixed 
points without bifurcation. In summary, combining the analytic discussion of section II with 
the numerical simulation we deduce that these maps are ergodic in certain region of their 
parameters space as explained above and in the complementary region of the parameters 
space they have only a single period one attractive fixed point, such that in contrary to 
the most of usual one-dimensional one-parameter or many-parameters family of maps they 
have only a bifurcation from a period one attractive fixed point to chaotic state or vice- versa. 

5 Conclusion 

We have given hierarchy of exactly solvable many-parameter family of one- dimensional 
chaotic maps with an invariant measure, that is measurable dynamical system with an in- 
teresting property of being either chaotic (proper to say ergodic ) or having stable period 
one fixed point and they bifurcate from a stable single periodic state to chaotic one and 
vice-versa without having usual period doubling or period-n-tupling scenario. 
Again this interesting property is due to existence of invariant measure for a region of the 
parameters space of these maps. Hence, to approve this conjecture, it would be interesting 
to find other measurable one parameter maps, specially higher dimensional maps, which is 
under investigation. 
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Figures Captions 

Fig.l. Plot of ^2^ a2 ( x ) ma P f° r a i = 1-5, a 2 = 2.9 ( curve with red color ) and plot of 
< ^2,2 Q2 ° < ^2, 1 2 Q2 ( a; ) ma P f° r a i — 1-5, «2 = 2.9, ( curve with green color ). The values of 
the maxima and minima (0 or 1) are independent from the parameters a±, a 2 , as shown 
in the figure. 

Fig. 2. Solid red surface shows the variation of KS-entropy Q^ 012 (x) , i n terms of the 
parameters a.\ and a 2 , dotted green surface curve shows variation of Lyapunov charac- 
teristic exponent of $2,2' a2 ( z ) > in terms of the parameters a\ and a 2 . 
Fig. 3. Solid red surface shows the variation of KS-entropy & 2 ^ a2 (x), in terms of the 
parameters a\ and a 2 , dotted green surface curve shows variation of Lyapunov charac- 
teristic exponent of & 2 ^ a2 (x) , in terms of the parameters a.\ and a 2 . 
Fig. 4. Solid red surface shows the variation of KS-entropy &£i a2 {x), in terms of the 
parameters a.\ and a 2 , dotted green surface curve shows variation of Lyapunov charac- 
teristic exponent of ^ 2 ,a2 (x), in terms of the parameters ol\ and a 2 . 
Fig. 5. Solid red surface shows the variation of KS-entropy 3>3 3' a2 (x), in terms of the 
parameters cti and a 2 , dotted green surface curve shows variation of Lyapunov charac- 
teristic exponent of $3 j 1 3 ' a2 (x), in terms of the parameters ol\ and a 2 . 

Fig.6. Bifurcation diagram of ^ 2 ] 2 a2 {x): for a\ > 0,a 2 > and a\ x a 2 < 4, it is 
ergodic and for a>\ x a 2 > 4, it has x — 1 fixed point. 

Fig. 7. Bifurcation diagram of $ 2 ^ a2 (x) and $ 2 ^ a2 (x): for a.\ > 0,a 2 > and 
«i x a 2 < 6, it is ergodic and for ot\ x a 2 > 6, it has x — 1 fixed point. 
Fig. 8. Bifurcation diagram of ^^ a2 (x): for a\ > x a 2 > | and a\ x a 2 < 9, it is 
ergodic, while for a± x a 2 < |, it has x = fixed point and a.\ x a 2 > 9, it has x = 1 
fixed point. 
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